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Abstract. - In the spinless Holstein model at half-filling the coupling of electrons to phonons 
is responsible for a phase transition from a metallic state at small coupling to a Peierls distorted 
insulated state when the electron-phonon coupling exceeds a critical value. For the adiabatic 
case of small phonon frequencies, the transition is accompanied by a phonon softening at 
the Brillouin zone boundary whereas a hardening of the phonon mode occurs in the anti- 
adiabatic case. The phonon dynamics studied in this letter do not only reveal the expected 
renormalization of the phonon modes but also show remarkable additional contributions due 
to electronic particle-hole excitations. 



Introduction. - The electron-phonon (EP) interaction leads in many low-dimensional ma- 
terials like MX chains, conjugated polymers, or organic transfer complexes [1,2] to structural 
distortions. Thus, Peierls transition and charge-density wave instability have been observed 
in such systems. Many interesting questions arise not only with respect to the associated 
metal to insulator transition but also concerning the signatures of the single-particle excita- 
tions in the different phases. The interest in models of electrons interacting with phonons has 
been renewed by intriguing findings which implicate an important role of the EP coupling in a 
wide range of materials with strong electronic correlations; high-temperature superconductors, 
manganites, or Cgo based compounds are well-known examples [3]. 

The one-dimensional Holstein model of spinless fermions (HM), 

H = -t£)(cta + h.c.) + woX! b t b i + 9^2(b\ + bi)ni, (1) 

is the simplest realization of a strongly coupled EP system and considers the local interaction 
g between the electron density n, = c\ci and dispersion- less phonons with frequency loq. The 
c\ (b\) are fermionic (bosonic) creation operators of electrons (phonons), and denotes 
the summation over all neighboring lattice sites i and j. With increasing g, the half-filled 
HM undergoes a quantum phase transition from a metallic state to a dimerized Peierls phase. 
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Because the HM is not exactly solvable a number of analytical and numerical methods have 
been applied (see references in [4]) to study the phase transition. In particular, the properties 
of the insulating phase have been shown to be sensitive to the relation between band-width 
and EP coupling [5]: a band insulator is found in the adiabatic case u>o/t <C 1 whereas a 
polaronic superlattice occurs in the anti-adiabatic limit ajo/t~^> 1. 

Recently, we applied the projector-based renormalization method (PRM) [6] to the HM 
at half-filling where both the metallic and the insulating case has been studied [4,7]. Here, 
we extend our work [4] to discuss the phonon as well as the electronic one-particle excitation 
spectrum in more detail. Most of the work is restricted to the adiabatic limit cuo/t <C 1 where a 
phonon softening is found if the EP coupling approaches the critical value g c of the transition. 
In contrast, a phonon stiffening occurs for uo/t 1. The most remarkable findings of our 
present work are the large contributions to the phonon spectrum which are caused by the 
coupling to electronic particle-hole excitations. Note that presently in the anti-adiabatic limit 
the insulating phase can not be studied because no stable solutions are found for t <C g. 

Theoretical approach. - We use the recently derived uniform description of metallic and 
insulating phases of the half-filled HM [4]. This approach employs the PRM [6] where an 
effective Hamiltonian TL — liniA^o Tlx is obtained by a sequence of unitary transformations, 
W(a-aa) = e Xx ' AX H\ e~ Xx - AX , by which transitions between eigenstates of the unperturbed 
part Tlo of the Hamiltonian are eliminated in steps. The respective transition energies are 
used as renormalization parameter A. Xa.aa has to be adjusted so that 7~t(\—A\) only contains 
excitations with energies smaller or equal (A — AA). In this way, an effectively free model, 

k>0,a k>0 9>0,7 

+ VNA b (b^ + h^+E, 

was found in Ref. [4] which is used here in order to calculate the phonon spectral function, 

1 f°° 

B(q,u) = ^ J 4>\]) e iut dt (3) 

and to consider the two electronic one-particle spectral functions 

K^) = ^ J" (<%(*) 4) ^ dt , A- k {u) = ^J" (4 c k (t)) e*« dt . (4) 

In Eq. J2J) a reduced Brillouin zone has been introduced in order to allow a dimerization of 
the system, and both the fermionic and bosonic one-particle operators as well as the model 
parameters have additional band indices, a,J3 = 0, 1. In Eq. @, <fi q = b q + b^_ q is proportional 
to the Fourier transformed lattice displacement. (ui) describes the creation of an electron 
with momentum k at time zero and its annihilation at time t whereas in (w) first an electron 
is annihilated. As is well-known, A^(u) [A^(w)] can be measured by [inverse] photoemission. 

To evaluate Eqs. © and Q within the PRM approach we use that expectation values are 
invariant with respect to an unitary transformation under the trace. Thus, B(q,cu), A^(ui), 
and A^ (to) can easily be computed if the phononic and electronic one-particle operators are 
transformed in the same way as the Hamiltonian (see Ref. [7] for more details). For instance, 
by taking the fully transformed phonon operator b q = <f q b^ q + fj q b- q + -4= ~J2 k ipk,qC k+q Ck we 



S. Sykora, A. Hubsch, and K. W. Becker: Phonon dynamics in the Holstein model 3 




-3 -2 -1 1 2 3 1 2 3 4 

ml t alt 



Fig. 1 - (Color online) Panel (a) shows the electronic one-particle spectral functions At(u>) (solid 
lines) and A^(ui) (dashed lines) for the metallic phase in the the adiabatic case, ujo = 0.4 1, g — 0.632 1 
where an one-dimensional HM with 500 lattice sites has been considered at half-filling. The Fermi 
energy is located at uj — 0. Panel (b) shows the corresponding one phonon spectral function B(q,uj). 



obtain for the phonon spectral function in the metallic region 

B{q,u) = %U(w - u t ) + ^-5{u + u>- q ) (5) 
+ 77 V m, q \ ; 0{e k+q -e h -u>) 

JV '—f Ek+q — Sk 

where terms with two bosonic creation or annihilation operators have been neglected. The 
coefficients <p q , f) q , and ipk,q can De determined within the PRM approach. Eqs. J3J and 
fulfill sum rules, f_ duj ui B(q,u>) = 1 and dcu [A^(cu) + A^(uj)] — 1, which are not 
affected by the PRM approach. All actual calculations are performed as described in Ref. [4] , 
and the half-filled HM is only considered in the one-dimensional case. 

Adiabatic Limit. - In the following we discuss the results for the one-phonon and one- 
electron spectral functions, B(q,uj), A^(uj), and A~^(uS). We start with the adiabatic limit 
(ivo/t <C I) where the metallic as well as the insulating phase can be described within the 
PRM approach. In panel (a) of Fig. ^ ( w ) anc ^ ^fc ( w ) are snown f° r different values of 
the wave vector k as functions of uj. The used parameters, lu = 0.4 1, g — 0.632 1, belong to 
the metallic region in the adiabatic limit so that no gap occurs at the Fermi level at u> = 0. 
Both functions are dominated by the coherent electronic excitation at ij- [8]. We also find 
incoherent excitations distributed over an energy range of about 2cj - 

As one can see in panel (b) of Fig.^ the frequency behavior of the phonon spectral function 
B(q,w) is also dominated by the coherent excitation at uj q [8] in this parameter regime. If 
the EP coupling g is close to the critical value g c of the metal-insulator transition, a strong 
softening of the coherent excitation at & q can be observed in the adiabatic limit for q values 
approaching the Brillouin zone boundary at q = tt. This phonon softening is caused by the 
coupling to electrons and has already been discussed in Refs. [4,7]. It will turn out that a 
phonon hardening occurs in the anti-adiabatic limit. 

Panel (b) of Fig. Q]shows that the energies of the incoherent excitations are distributed over 
an energy range of the order of the electronic bandwidth. According to Eq. (J5J, the energy 
spreading is caused by the energy differences of electronic particle-hole excitations (sk+q — £k) 
where k runs over the whole Brillouin zone. From the numerator in Eq. [/(sfc) — f(ik+ q )], 
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Fig. 2 - (Color online) As Fig. 0but for u>o = 0.1 1, g — 0.3 1 (adiabatic case, insulating phase). 



one concludes that either ik is smaller and ik+q larger than the Fermi energy or vice versa. In 
agreement with Fig.^ Eq. JHJl also shows that the energy range of the incoherent excitations 
increases with increasing q. Due to the absence of an electronic gap in the metallic phase, 
incoherent excitations with small energies are present for q- values close to q — and q = tt. 

Next, let us consider the ordered insulating phase in the adiabatic limit, u>o/t <C 1, where 
ujq = 0.1 1 and g = 0.34 t have been chosen. The gap of order A w 0.37 1 in the electronic 
excitation spectrum [compare panel [a] of Fig. Ej clearly shows that the system is in the 
insulating state. Again the electronic spectral functions A^{uj) and A~^(ui) are dominated by 
the coherent excitation at ik, and the spreading of the incoherent excitations is similar to that 
of the metallic state [compare with panel (a) of Fig. As one can see in panel (b) of Fig. |5J 
the spectral weight of the coherent excitation is for the phonon spectral function B(q, u) again 
huge in comparison to that of the incoherent excitations. However, it is important to notice 
a significant difference to the metallic state [compare panel (b) of Fig. No incoherent 
excitations occur at low energies which corresponds to the gap in the electronic spectrum. 

Anti- adiabatic Limit. - Next we discuss the anti-adiabatic limit where u>o/t 3> 1 holds. 
Unfortunately, in this limit we are restricted to the metallic case g < g c , because no stable 
solution for the insulating phase g > g c is found within the PRM approach. The reason is 
that the absolute value of the EP coupling g is too large in this case. As one can see in 
Rcf. [4], the renormalization equations were derived by starting from an uncorrelated model 
by successively eliminating high-energy excitations which prevents the application of the PRM 
scheme for extremely high coupling parameters g. The electronic spectral functions At(ui) 
and A~j~(u>) are shown in panel (a) of Fig. for the metallic phase, where ujq = At and 
g = 2.82 1 have been chosen. The missing gap in the electronic spectrum reveals the metallic 
state of the system, and the k dependence of the sharp coherent excitation peak corresponds 
to the dispersion ik of the electrons. The coherent excitation is separated from the incoherent 
excitations by an energy of the order of the phonon energy lu . Here, a renormalized one- 
electron creation operator ct(A) has been used which includes phonon operators only in linear 
order (see Ref. [7] for details). This coupling leads to the energetic separation of order u)q. If 
operator terms with more than a single phonon operator were taken into account additional 
excitations at higher frequencies (2woj 3wo, etc.) would appear. This is confirmed by recent 
ED results of [9]. 

The metal-insulator transition in the anti-adiabatic limit can be understood as the for- 
mation of small immobile polarons which arc electrons surrounded by clouds of phonon ex- 
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Fig. 3 - (Color online) As Fig. 0but for u>o = 4t, g = 2.82 i (anti-adiabatic case, metallic phase) 



citations. In the PRM approach polarons correspond with the fully renormalized electronic 
excitations ik- If one compares the renormalized electronic band width for the adiabatic case, 
which can be read off from the coherent excitation in panel (a) of Fig. ^ with that of Fig. 
in the anti-adiabatic case one observes a strong reduction of the band width. This reduction 
clearly indicates localization tendencies in the system which might allow the determination of 
the critical EP coupling g c of the metal-insulator transition in the anti-adiabatic limit within 
our PRM approach [10]. 

Panel (b) of Fig.[3|shows the phononic spectral function B(q, u>) in the anti-adiabatic limit. 
It is important to notice three significant differences to the adiabatic case [compare with panel 
(b) of Fig. Firstly, the coherent phonon excitation no longer shows a softening behavior, 
instead an hardening of the phonon modes is observed from a value ojq at q = to higher 
energies for all q > 0. Secondly, the incoherent contributions from particle-hole excitations 
have gained considerable weight as compared to the adiabatic case [see panel (b) of Fig. 
The energy range of the incoherent excitations is again of the order of the electronic band 
width. Finally, a huge elastic contribution is found in B(q,u>) for q = tt. 

Discussion. - Phononic and electronic spectral functions for the one-dimensional HM 
at half-filling were recently evaluated in Ref. [9]. In this work, the same parameter values 
as used here were studied by exact diagonalization techniques (ED) for system sizes of up 
to 10 lattice sites. The results for the electronic spectral functions (ui) and A^(uj) agree 
quite well with our analytical results in the whole metallic regime. In the phonon spectrum 
two essential differences are found. In the ED-approach [9] a strong excitation is found at 
frequency oj — for wave vector q — both in the adiabatic and in the anti-adiabatic case. 
In contrast, the results in Figs. ^and[3]show excitations for small q only for finite frequencies 
oj > 0. This additional zero- frequency peak in Ref. [9] is caused by a different definition of 
the phonon spectral function. Instead of the commutator spectral function B(q,u>), in [9] a 
phononic correlation function, B{q,uj) = f-oo (QqW) ^J) e %ut dt, is investigated so that an 
additional zero-frequency contribution from the center of gravity motion contributes, which 
is dropped in the commutator spectral function as defined in Eq. 

The second difference between our results and those of the ED approach of Ref. [9] is more 
important: In the phonon spectrum of [9] almost no small frequency excitations are found in 
the vicinity of q — tt whereas the present PRM approach for B(q — tt, uj) in the anti-adiabatic 
limit gives a strong zero-energy peak at q — tt. For a more detailed discussion let us consider 
the phonon spectrum for different q- values in the vicinity of q = tt [panel (b) of Fig. [!]• 
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Fig. 4 - Low energy range of the phonon spectral function as shown in Fig. [3] for small q values 
[panel (a)] and for q values close to the Brillouin-zone boundary [panel (b)]. 



As follows from (JjJJ, the incoherent excitations in B(q,u>) are caused by the coupling of the 
phonons to electronic particle- hole excitations with energies \£k+q — £/c| where the wave vector 
k runs over the whole Brillouin zone. The energies £k+q and ik have to be below and above 
the Fermi level with momentum kp. Thus, only electronic particle- hole excitations contribute 
with either |fc| < kp and \k + q\ > fcp, or |fc| > kp and \k + q\ < kp. For q = ir this condition 
is fulfilled for all fc-values of the Brillouin zone, and incoherent excitations are found in an 
energy regime of the order of the electronic band width. For q = ir the smallest possible energy 
difference \£k+q — £k\ is the same as the energy difference of two neighboring k points around 
kp with the smallest possible distance in fc-space. This distance is of order 1/N and might be 
not small for a small system size. This explains the absence of low-energy excitations in the 
phonon spectrum of Ref. [9] close to the Brillouin-zone boundary (where only system sizes 
of order N=10 have been considered). We could confirm this mechanism within our PRM 
approach: If we reduce the system size from 500 to 66 lattice sites the low-frequency peak in 
B(q = it,uj) is almost completely suppressed, and its position is shifted to higher energies. 

Let us also discuss the case of small q- values q <C 7r, where only k values from the sum 
in (jSJ) can contribute which are located in a small wave vector region around the Fermi 
momentum kp. Assuming a linear momentum dependence for energies close to Fermi energy, 
the incoherent electronic excitation energies can be replaced by \ik+q — £k\ ~ vpq. [yp is 
the Fermi velocity.) For small q the position of the incoherent excitations is proportional to 
q, as seen in panel (a) of Fig. On the same time the intensity of the incoherent peak is 
approximately q independent which follows again from Q: The linear q dependence of the 
denominator is canceled by the number of contributing k points which is also proportional to 
q so that an approximately q independent intensity results. 

Finally, we use the PRM to evaluate the Luttingcr liquid parameters u p , K p in the metallic 
phase at half-filling. It is commonly accepted that the spinless one-dimensional HM belongs 
to the Tomonaga-Luttinger universality class, and at T — a Luttinger liquid should be 
realized. Due to the quite different behavior of the phonon spectral function in the adia- 
batic and in the anti-adiabatic regime one expects different results in both cases. For the 
velocity u p of the charge excitations, we use a finite-size scaling of the energy gap [11,12], 
A p = ^(Eq 1 + Eq 1 — 2E ) = -rf , where N is the lattice size, and E denotes the ground- 
state energy of the half- filled system. E^ 1 is the ground-state energy of the system with 
±1 fermions away from half-filling. To compute the effective coupling constant K p we take 
the first derivative of the charge correlation function C{q) at q = 0, dC Q^ | g =o = where 
C{q) — jf X)i(=>%- e lq ^ Ri ~ R ^{SniSnj). Results are given in Table [I] In the anti-adiabatic 
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Table I - LL parameters in the adiabatic and in the anti-adiabatic regime. To account for the missing 
spin degree of freedom, K p has been multiplied by a factor of 2. 



regime the parameter u p for the kinetic energy is strongly reduced when g is increased which 
corresponds to the band width reduction as discussed above [compare panel (a) of Fig. 
This behavior of u p is very similar to recently found results [13]. In contrast, the results for 
K p represent an unsolved puzzle: Here, K p is always smaller than 1, whereas in Ref. [13] K p 
was smaller than 1 only in the anti-adiabatic regime but larger than 1 in the adiabatic regime. 

Summary. In this paper we have considered the phononic and electronic spectral 
functions of the one-dimensional HM at filling. In particular, we have studied the spectral 
signatures of the metallic phase both in the adiabatic and the anti-adiabatic limit, and of the 
insulating phase in the adiabatic limit. (The insulating state could not be investigated in the 
anti-adiabatic limit because the employed PRM approach breaks down in this case.) We find 
a quite strong coupling of the phonon dynamics to electronic particle-hole excitation which 
leads to incoherent contributions to the phononic spectral function. In particular, a dominant 
incoherent low-energy peak is observed for the metallic phase in the anti-adiabatic limit. 
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